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SUMMARY 


An equivalence between temperature and pressure, heat conduction and 'stiffness', 
and heat capacity and 'mass' is defined, enabling acoustical modal analysis of 
arbitrary three-dimensional volumes. The transient heat transfer analysis rigid 
format in NASTRAN, RF9, has been ALTERed providing the acoustical analysis capa- 
bility. Examples and ALTERs are included. 


INTRODUCTION 


A twenty-node isoparametric acoustic finite element model was developed in 
Reference 1 for analyzing the acoustic modr of irregular shaped cavities. In the 
present paper, 

1. the existence of an identical element in NASTRAN (IHEX2) and, 

2. the recognition of the similarity between the acoustical matrices of 
Reference 1 and the thermal matrices of NASTRAN 

have enabled the posing and solution of the acoustics eigenvalue problem for arbi- 
trary three-dimensional cavities bounded by 'hard' acoustic surfaces. A simple 
modification of the transient heat transfer rigid format RF9 provides the acoustics 
analysis formulation in NASTRAN. 

SIMILARITY OF ACOUSTICAL AND THERMAL MATRICES 


Summarizing the finite element formulation of Reference 1, the pressure p in a 
volume V bounded by a surface S satisfies the three-dimensional wave equation and 
boundary conditions: 


V 2 p + (w 2 /a 2 )p 
o 


Vp • n 


0 in V, and 

(1) 

0 on S, 

(2) 


where uj is the frequency of vibration of the acoustical mode, ao is the speed of 
sound, and n is the outward normal to S. Representing the volume V by an assemblage 
of three-dimensional finite elements, the corresponding eigenvalue problem becomes 

[K - io 2 M] {p> - 0 , (3) 

where for the ith element 

[k] t - / CBl^ [B] ± dV , (4) 

V i 
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(5) 


[mJ i “ / ^i- N -£ LN J i dV ’ 


3 

"3x 


[B] 


3 _ 

3y 



3 _ 

3z 


( 6 ) 


The shape function LNj^ approximates the pressure within the ith element in terms of 
the nodal pressures {pl^ as 

p = Ln^ {p> i . (7) 

A comparison with the NASTRAN heat transfer analysis capability (Ref. 2) 
indicates that, 

1. equation (4) is identical to the heat conduction matrix, if the thermal 
conductivity is unity, and 

2. equation (5) is identical to the heat capacity matrix, if the thermal 
capacity per unit volume is 1/ao. 

The temperature degrees of freedom in the thermal analysis are taken to corre- 
spond to the pressure degrees of freedom in the acoustics analysis. 

The above correspondence together with the real eigenvalue analysis module 
'READ) permit the determination of the acoustical modes and frequencies. An ALTER 
package to be used in HEAT APPROACH RF9 is included in the Appendix. 


EXAMPLES 


The two examples of Reference 1, shown in Figures l and 2, were analyzed using 
the modified RF9. The results, generally in agreement, are presented in Tables 1 and 
2, and Figure 3. It is noted that the missing geometric dimensions in Figure 2 were 
scaled from the figure of Reference 1. 


CONCLUDING REMARKS 


With a simple modification, and a 'redefinition' of thermal conductivity and 
capacity, the transient heat transfer rigid format in NASTRAN has been used to deter- 
mine the acoustical modes and frequencies of arbitrary volumes. The volume can be 
modelled by any of the solid elements permitted by RF9. Although only acoustically 
hard surfaces have been considered in this paper, simple extensions to other boundary 
conditions are considered to be possible. 
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* ALTER HEAT K'F 9 (NASTRAN RELEASE APRIL 1982) 

* 

ALTER 26 1 26 $ SET N0MGG- 1 BECAUSE C r ERROR IN I HEX 2 LOGIC ♦ 

EMG WEST » CSTM r MPT r HIT » GE0M2 , /HKELM r HKDICT r DUM1 1 DU M2 r HBELM r HBDICT/ 

SfNrNOKGGX/ 1/ SrNrNOBGG $ 

ALTER 67,67 $ COMPUTE AND PRINT MODES 

DPD DYNAMICS > GPL , HSIL > HUSET/GPLD r HSILD , HUSETD * TFPOOL » HDLT r » » 

HNLFT > HTRL r FED > H E 0 D YN/HLUSET/S r N - Hi. USE I'D/ 1 23/S , N » NODE I / 

123/1 23/S , N ? N0NI..F T/S?N» NOTRL/S » N t NOEED//S > N , NOUE $ 

PARAM //*MPY*/NEIGV/1/-1 * 

READ HKAA > HBAA * > , EED r HUSET » CASECC/LAMA * PHI A, , OEIGS/ 

*MODES*/S»NfNEIGV * 

OFP OEIGS r LAMA » , , , // $ 


SDR. I 
SDR2 

OFP 

PLOT 


PRTMSG 


HUSET , , PH I A , , r!iuO»GM» , , , /PHIG , »/t/#REIG# ♦ 

CASECC , CSTM , MPT , D I T » HEQDYN » HSIL » » » BGPDP , LAMA r rPHIGr 


HEST , , / , , OPP IG , 
OPH .I.G, f f „ , „ $ 


»»PPHIG/#RETG* 


* 


PLTPAR r GPSETS » EL SETS , CASECC , BGPDT t HEQEXINr HSIP , , PPHIG , HGPECT , 
/PL0TX3/HNSII./HLUSEP/JUMPPL0T/PLTFLG/PFILE * 

PL0TX3// $ 


ENDALTFR 


; 


TABLE 1. 

ACOUSTIC FREQUENCIES OF A RIGHT-ANGLED 

PARAI.LELOPIPED 


(a) 

Symmetric Modes 


Mode 

Exact Frequency* 

Ref. 1 

This Paper 

o 

o 

699 

702 

680 


l 

to 

O 

o 

1398 

1542 

1495 

o, 0, 1 

1500 

1506 

1460 

* 1, 0, 1 

1655 

1664 

1613 

2, 0, 1 

2057 

2155 

2090 

3, 0, 0 

2097 

2525 

2448 

3, 0, 1 

2579 

— 

2925 
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TABLE 1. ACOUSTIC FREQUENCIES OF A RIGHT-ANGLED PARALI.ELOPIPED (ryc/T) 
(Contd. ) 

(b) Antisyrunetric Modes 


Mode 

Exact Frequency* 

Ref. 1 

This Paper 

0 , 1 , 0 

1250 

1255 

1255 

1, 1, o 

1432 

144C 

1430 

2, 1, 0 

1876 

1988 

1952 

0, 1, i 

1953 

1963 

1928 

1, 1, 1 

2074 

2089 

2048 

2, 1, 1 

2401 

2499 

2440 

3, 1, 0 

2442 

2875 

2808 

3, 1, l 

2866 

3346 

3233 

* £ l. 

. iLSL i(— ) 2 + 
m, n 2 U «, ' 

(f- )2 + » 

y z 

(Ref. 1) 

V 

l , l sides of the 

y z 

parallelepiped 



TABLE 2. 

ACOUSTIC FREQUENCIES OF MODEL VAN 

(cyc/T) 


(a) Symmetric 

Modes 



Ref. 1 



Mode 

Experimental 

Calculated 

This Paper 

1 

606 

593 

562 

2 

1174 

1150 

1104 

3 

1549 

1556 

1 524 

4 

1613 

1605 

158(/ 

5 

1817 

1829 

1751 

5 

1992 

2026 

1993 



(b) Antisymmetric 

Modes 


Mode 

Pef . 1 

Experimental 

Calculated 

This Paper 

1 

1220 

1168 

1118 

2 

135/ 

1317 

1245 

3 

1673 

1634 

] 568 

4 

1996 

} 956 

1940 

5 

2021 

1 997 

1968 

6 

2176 

2187 

21 09 


* 
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FIGURE 1. DETAILS OF RIGHT-ANGLED PARALLELOPIPED 
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